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Abstract 

The electrostatics properties of composite materials with fractal geometry 
are studied in the framework of fractional calculus. An electric field in a com- 
posite dielectric with a fractal charge distribution is obtained in the spherical 
symmetry case. The method is based on the splitting of a composite volume into 
a fractal volume ~ r d with the fractal dimension d and a complementary host 
volume Vh = V3 — Vd- Integrations over these fractal volumes correspond to the 
convolution integrals that eventually lead to the employment of the fractional 
integro-differentiation. 

Keywords: Fractional integral, Fractional derivative, Fractal geometry, 
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1. Introduction 

The electrostatics of composite materials is a ubiquitous and fundamen- 
tal phenomenon with a variety of applications in physical, electro-chemical, 
geophysical, and biological systems, see e.g. The investigation of the 

electric properties of composite materials has been a long-lasting task, start- 
ing, e.g., from a calculation of the permittivity of a dispersed mixture [3| and 
continuing with the modern investigations of nanosystems, for example, of an 
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enhanced scattering of electric fields in nanosystems |4|, Q] . The concept of 
fractals and fractal geometry is widely employed in these studies, and leads to 
essential progress in our understanding of the physical properties of compos- 
ites. A wide class of the physical applications of fractal structures has been 
reviewed well, see e.g., In its turn, the fractal concept makes it 

possible to involve the application of the powerful methods of fractional integro- 
differentiation. An illuminating example of this relation is fractional diffusion 
on a comb [t], floL [ill ] , where the relation between the comb geometry and non- 
Markovian transport is straightforwardly established in the fractional Fokker- 
Planck equation [l(J [ill Hi ■ As obtained in [3| (see also [l] chapter 39), the 
mathematical formulation of the transform of the Brown line-to-line function 
B(t) (continuous, but not differentiable) into fractional diffusion, described by 
fractional line-todine function B a (t) with < a < 2, is expressed by means of 
the fractional integration [j] 

B a (t) = -±- f {t-sy- l dB{s). (1) 

The concept of the integration and differentiation of non-integer orders already 

nn 

arises from the works of Leibniz, Liouville, and Riemann, see e.g., 

\m,\m- Its 

application is related to random processes with power law distributions. This 
corresponds to the absence of characteristic average values for processes that 
exhibit many scales [l?! [3] . 

By analogy with Eq. ([1}, a link between fractal geometry and fractional 
integro-differentiation [19J is constituted in the procedure of averaging an ex- 
tensive physical value that is expressed by means of a smooth function over a 
Cantor set, which leads to fractional integration. In its eventual form, it has 
been presented in Ref. 20j. The main idea of filtering or embedding a matter 



inside a fractal is the construction of a convolution integral 

F(x) = W(x)*f(x) = [ W(x-y)f(y)dy, (2) 



1 Note that fractional Brown motion was originally introduced by Kolmogorov in 1141 and 
was later rediscovered and explored in greater detail by Mandelbrot and Van Ness in 
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where the function W(x) obeys the scaling relation W(x) — ^W{bx) with a 
solution W{x) — x a A(x). Here a = y^r and A(x) = A(bx) is a log-periodic 
function with a period In b. When this scaling corresponds to a Cantor set with 
a = 2, and A(x) is defined explicitly, one performs averaging over this period 
|2o| | and obtains Eq. ([2]) in the form of a fractional integral 

(F{x)) = (x*- 1 * /(*)) = V(a) IZf(x) , (3) 

where V(a) — 2 h ^ 2 X and T(a) is a gamma function, while qI" designates the 
fractional integral of the order of a in the range from to x: 

= ^t/V- xT-'fWx' . (4) 
r (") Jo 

Analogously one can define the fractional derivative as the inversion operation: 
D" — q/jT". Extended reviews on fractional calculus can be found in [15I [if]. 



211 ] . for example. See Appendix A. 

This mathematical construction is relevant to the study of electrostatics of 
real composite structures. The main objective of the present research is to 
solve a standard electrostatic problem, namely, the derivation of the electric 
field in a composite dielectric with a fractal charge distribution. This problem 
was first considered in 22J, in a slightly different mathematical formulation, 



where integration was carried out over the fractal with fractal dimension a and 
infinitesimal fractal volume dr a inside a homogeneous dielectric. 



2. Fractional Distribution of Charges 

First, we demonstrate this fractal property of embedding a matter for elec- 
trostatics by implementing it in Maxwell's integral equation for the electric field 
E. We formulate the problem in terms of a characteristic function. Let any 
charge distribution on a fractal be F a (r), where r = (x, y, z) is described by the 
characteristic function x( r )i such that x(r) = 1 for r £ F a (r) and x(r) =0 for 
r ^ F a (r). Therefore, the fractal filtering of Eq. ([3]) now reads 

X (x')f(x')dx' -Lr [ X (x - x'y^fWx' = IZf(x) . (5) 
r (") Jo 
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We obtain the same integration for the y and z coordinates. 
The Gauss theorem in Cartesian's coordinates reads 

j V • E dxdydz = olWi V ' E = Q • ( 6 ) 

where an electric charge Q with a fractional distribution is embedded in a three- 
dimensional volume. Without restricting the generality, we consider the fractal 
charge in the first octant of the Euclidian space. Therefore, we have from Eq. 
(5) 

Q = Q(x, y, z) = IxX( x )ollx(y)ollx( z )p( x , V, z ) => oI"oIyoIzP(x, y, z) , (7) 
where p(x, y, z) is the charge density, and the constant multiplier V(a) is inserted 



in p; for an alternate definition see also |22| . The charge in Eq. ([7]) corresponds 
to a fractal structure that is the direct product of the Cantor sets defined in Eq. 
©, namely (F a {x)) x {F a (y)) x (F a (z)). 

To proceed, we obtain the correct dimensions of both the charge and the 
electric field, respectively. To this end it is convenient to work in the dimen- 
sionless space variables — » x, — y y, — > z, and the properly defined 
fractal charge reads from Eq. ([7]) 

Q(x, y, z) = I2 I« I?l 3 p(x, y, z) . (8) 

The conventional integral in Eq. ([B]) yields the dimensionless integration qI^qI^qI^ 
loolxolyoll an d dimensionless differentiation V — > ^V. Here Iq is the charac- 
teristic size of a physical fractal, which is always finite. For example, it can 
be the minimal scale of self-similarity of a physical fractal. In what follows we 
always bear in mind that the space variables are dimensionless, and the charge 
density is properly scaled such that fractal charges and electric fields have cor- 
rect dimensions. 

When the charge density is a constant function, the fractional convolution in- 
tegrals become a simple integration over a dimensionless volume (xyz) a /T 3 (a + 
1). Comparison of Eqs. © and (UJ) yields Maxwell's equation in dimensionless 
space coordinates 

d*Q (xyz)*- 1 
V ' E - dxdydz --THaT loP - (9) 
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It has a solution 

E = -rtt [*a a (y*) a_1 + WW*' 1 + iz^xyY- 1 ] . (10) 

al J (a) L J 
When the charge density p(x,y,z) is not a constant, the integrands in Eqs. ([6]) 
and ([7|) cannot be straightforwardly compared to each other. Therefore, we ap- 
ply the fractional differential qD® to both sides of Eq. ^ . Taking into account 
that D^ I°f(x) = f(x) and D% I%f(x) = D^f(x) (see Appendix, Eq. 
(|A. 6p ). one obtains Maxwell's equation in form either 



/ x 1 - < V y 1 - ( Vi- a V • E(z, y, z) = p(x, y, z) , (11) 

or 

V • E(.t, y, z) = a Dl- a Dl- a Dj l - a p{x, y, z) . (12) 
For a constant p it reduces to Eq. ([9]). 



3. Example: Electric Field of an Infinite Fractal String 

Now we consider an example where the filtering by means of the convolution 
integral in Eqs. © and Q can be considered as a superposition of electric 
fields, that is the fractal filtering of the electric field. Let us find an electric field 
of a charged fractal string with the fractal dimension < v < 1 and a constant 
linear charge density p\ . We use cylindrical symmetry and describe the problem 
in terms of the string coordinate z and the distance from the string r. The E r 
component of the electric field at point (r, z) of an infinitesimally small charge 
with coordinate z' is dE r (r. z) = prrdz where we used that the linear 

io[r 3 +(2-*') 2 ] 3 

charge density is scaled ^ . Therefore, the complete electric field at this point is 
a superposition of the charges along the fractal string. This yields the fractional 
integration 



/ 1» IJ-oc [r 2 + (z - z') 2 ]i Jz [r 2 + (z - z'Y]i J 
Performing the variable change y = (z — z')/r, one obtains for the electric field 

E r (r,z)^E r (r) = 2 -^, (14) 
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where A(v) = ^tt-t f_°° d 4-. When v = 1, this yields a well known result 

for the electric held of an inhnitely long line charge. In dimension variables 
the electric held reads E r {r) = 2A{y)(^\ i.e., the electric field of a 

fractal string decreases with the distance r (r/Zo) 1_l/ times faster than one of 
the continuous counterparts. 



4. A Case of Spherical Symmetry 

Now we obtain an expression for the fractal charge in spherical coordinates. 
We define a charge on a random fractal, which has a homogeneous distribution 
in a three-dimensional medium with charge density po- A fractal mass inside 
a ball of radius r is A4(r) ~ r d . As shown in Sec. 2, the convolution in Eq. 
(J7J) determines the map of the charge into the fractal volume. In Cartesian's 
coordinates it corresponds to the fractional Riemann-Liouville integral with el- 
ementary fractional volume 21, 3| 

I xy z I s ^ 

dVd ^YHd/W dxdydz - 

In the spherical coordinates, which correspond to the Reisz definition of the 
fractional integral, the elementary fractional volume is 

dV d = — ^ ; |r| rf ~ 3 r 2 dr sin OdOdcj) . 

In the spherically symmetrical case, the convolution integral (0 takes a simple 
form. Hence, the random fractal can be considered as a direct product of the 
two-dimensional sphere 5*2 and a Cantor set F a of the dimension a, such that 
the fractal dimension of S% <g> F a is d = 2 + a 24]. Therefore, the mapping of 



the charge inside the fractal is determined by the convolution 

Q{r) = ^-r f(r - r')^ 1 P o(r')r' 2 dr' , (15) 



r(«) Jo 

where all constants are taken inside the fractional charge density p. In what 
follows the charge density po is a constant value. This choice immediately yields 
an expression for the charge from Eq. (|15[) 



G 



that correlates with the fractional mass M(r). This also yields the average 
fractal charge density of the entire composite 

rfr)~^~*r--». (17) 

This expression yields the corect result for the electric field inside the fractal 
dielectric 

m = J&W^ poRd - 2 (18) 

when R < r [2^] . Note also that for the electric field outside the charged 
fractal, when R > r, we have for the average fractal charge density p{R) — > 
p(r')Q(r' - r). This yields for the electric field E(R) = Here Q(z) is the 

Heaviside step function. 

5. Electric Field in a Fractal Dielectric Composite 

The situation with respect to permittivity is different and needs special care. 
Let us consider the Maxwell equation (|6]) inside a volume that contains a fractal 
dielectric. We take into account that the electric field is in a random fractal of 
a volume Vd ~ r with permittivity s\ and in a free space of a complementary 
host volume Vh = V3 — Vd with permittivity e 2 . We suppose that the random 
fractal is a direct product of the two-dimensional sphere S 2 and a Cantor set 
F a of the dimension a, such that the fractal dimension of S 2 <8> F a is d = 2 + a. 
We also suppose the spherical symmetry, \7eEg — V^i?^ = 0. The electric 
displacement D(r) = e(r)E(r) is not differentiable, since the permittivity is a 
discontinuous function on the fractal: 

n / £U if r G V d 
e(r) = t (19) 

[ e 2 , if r £ V h . 

The permittivity of a two-phase composite in Eq. (1191) can be defined by the 
characteristic function e(r) — £ix(r) + £2[1 — x( r )]- The Maxwell equation for 
the f component reads 

(ei - e 2 )x{rWrE r + e 2 V r E r + (ei - e 2 )E r \7 rX (r) = p . (20) 
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The last term in the l.h.s. corresponds to a charge, which is due to the polar- 
ization of the fractal and to the discontinuity of the electric field on the fractal 
interface [3[. Now we use the Gauss theorem, which for the chosen geometry 
reads 

Q(r) = 4tt f \7 r D r {r')r' 2 dr' = 47r / r V 2 V r D r (r) . (21) 
Jo 

Taking into account that oIrX(r)r 2 f(r) ~~ ^ of? f 2 f{r), we obtain from Eqs. (|20|) 
and (ED 



-^Q(r) = {e 1 -e 2 ) I^r 2 V r E r {r)+e 20 iyV r E r {r) + {e 1 -e 2 ) a I 1 r E r {r)r 2 V r x{r) . 

47T 

(22) 

The last term in Eq. (|22|) can be written in the form (see Appendix B) 
oI^E r (r)r 2 . Finally, taking into account Eq. p^|) , we obtain the integral 
Maxwell's equation for the electric field in the form 

{e l -e 2 ) iy 2 ^ r E r {r)+e 20 iyV r E r {r)+ {ei ~ £2) \ l^E r {r)r 2 = -^— r 2+a 

£ 2 1 {a + 3) 

(23) 

We take into account that r 2 \7 r E r (r) = d r [r 2 E r (r)] = d r G(r), where we in- 
troduce a new function, G{r) = r 2 E r (r). We also introduce a dimensionless 
parameter 

e=^^. (24) 
£2 

Applying the fractional derivative o-D" to both sides of Eq. (|23|) and taking into 
account that ^E>%r v = p^^z^ r q ~ d one obtains 

ed r G(r) + Q ll- a d r G{r) + e 2 G(r) = ^ . (25) 

£2 

Here 

l}- a d r G(r) = f(r - r')- a d r -G(r')dr' . (26) 

r U - a ) Jo 

It is worth noting that for a homogeneous medium with a fractal distribution of 
charges, when Si — e 2 , Eq. (1251) reduces to Eq. ([9]) in the spherical symmetry. 

5.1. Solution by Expansion 

Eq. (|25|) is the linear fractional integro-differential equation that can be 
solved by the Laplace transform. Defining the Laplace transform C[G(r)] = 
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G(A), one obtains the solution of Eq. (J25]) in the form 

G(X) = ■ \ =r . (27) 

v ' e 2 A 3 [A + eA + e 2 ] y ' 

Expanding the denominator near eA, one obtains 

9/1 00 / i\n " ,2m 

(5(a) = ^ y y c™——^—, — - , (28) 

n— m— 

where C™ = r ( m + 1 jr (n- m + 1 ) ' Carrying out the inverse Laplace transform, one 
obtains for the electric field 

9„„. 00 ( 1 \n " 2m n(l—a)+ma 

Mr) = ^ y LiL y Cr rr ; r . — 1 . (29) 
£?e e 1 n l-a +Hitt + 4 

n— m— J 

When e< 1, one takes only m — in the second summation. This yields 

n— J 

The sum in Eq. ([30| is the definition of the Mittag-Leffler function (see e.g., 

HQ) 

°° r_ (l-a)/ ]n 

%- Q , 4) (-^/e)=y r ; n(1 _ a) / ; 4] . (3D 

In the limit e — >• 0, the asymptotic behavior of the Mittag-Leffler function for 
^ 3> 1 has a power law decay (see e.g., H^H) 

W-^-ror^. (32) 

and we arrive at the solution (TTH]) (see also 22[) 

^ = i^T3)^- < M > 

Since the charge is contained inside the fractal, we also used here that e > 0. 
This condition also ensures the transition to a homogeneous dielectric with 
£1 = £2 and is valid for a — 1, as well. 

In the opposite case when r< - - <C 1, the solution for the electric field reads 
EJr) ~ ^mi. 
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5.2. Approximations for e< 1 and e> 1 

Note that the expansion (|28p in not always valid. Therefore, it is convenient 
to use an approximation inside the denominator in Eq. (1271) . In the limiting 
cases t <1 and e 3> 1 it can be simplified. In these cases one can look for the 
solutions using an expression for the Laplace transform for the Mittag-Leffler 
function for two parameters £( a ^(—az a ). It reads [l^ 

e-^z^S ( „ )j8) ( Ta8 ") d* = y - - , (Re(i/) > |o| - ) . (34) 

We start by considering e <C 1. Neglecting e 2 , we have 

G(A) « 2 f° ■ p- 1 - . = (35) 
£ 2 eA 3+Q [A 1 -" + 1/e] 

Comparing Eqs. (jM)) and (|35p one obtains 

Br(r) « ^rf (1 _ a , 4) (-^) > (36) 



which yields the same expression for the electric field as in Eq. (|30|) . 

In the limit e 3> 1, we neglect A Q in the denominator in Eq. (|27[l. which 
yields the following solution for the electric field 

E r (r) w — r£(i 4)(— er) . (37) 

Note also that the dimensionless radius can be both r ^ 1/e and r <C 1/e. 
Therefore, for tr > 1 we use the asymptotic behavior of the Mittag-Leffler 
function (|3"2j> and obtain the asymptotic behavior of the electric field E r (r) rs 
^j*r-. In the opposite case, er <C 1, we have £(i- a ,4) ~ 1, which yields E r (r) rs 

gPo r 

5.5. Direct Product of Cantor Sets 

In the preceding section, the random fractal was considered as a direct prod- 
uct of a two dimensional sphere S2 and a random fractal f a of the dimension 
a, such that the fractal dimension of S2 <8> F a is d = 2 + a. It is worth no ting , 



that this is only one of many possible constructions of a random fractal 



24|. 



Therefore, it is instructive to consider a more general example. Let us consider a 
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fractal dielectric as a product of the random Cantor sets F a (x) x Fp(y) x F 7 (z), 
such that the electric charge due to Eq. ([7]) is 

Q{x,y,z) = —— — — — -, x,y,z.O 

r(i + a)r(i + p)r(i + 7) 

where a + f3 + 7 = d is the fractal dimension <i < 3. We simplify the analysis 
by neglecting the effect of polarization of the composite. This simplification is 
valid only in the case when e <IC 1. As shown, the polarization charge is of the 
order of e 2 , and it can be neglected in the Maxwell equation. Therefore, the 
integral Maxwell's equation reads now 

Q(x, y, z) = ( El - e 2 ) qI^o^V • E(x, y, z) + £ 2 o4o^ViV • E(x, y, z) . (38) 

Applying the fractional derivatives oD" D^ DJ to both sides of the equation 
and denoting F(x, y, z) — V • E(ir, y, z), one obtains the Abel equation 2l|, |26| 

^ = eF(x, y, z) + J^V^" Vi" 7 ^, 2/, z) ■ (39) 
£2 

Here -F (x, y, z) is an effective charge density. We solve this equation by a method 
of consequent approximations [2^. The first step is F{x, y, z) = j^+Fi(x, y, z), 
which yields the following equation for F\(x, y, z) 

= eF 1 {x,y,z) + F 1 (r) + o^o^o^ 7 ^, V,z) 

+ V - (40) 

e 2 er(2-a)r(2-/3)r(2- 7 ) 1 ; 

Then, defining F x {x, y, z) = - j^ r (2-I)r(2-ff)r(2- 7 ) + F *( x i V' z ) and repeating 
the same procedure for finding F n (x, y, z), we have the solution for the effective 
charge density in the form of the sum 

00 

F(x,y,z) = — ^—y2F n (x,y,z), (41) 

where 

\-x 1 - a y 1 -f i z 1 -i /(e)] n 



r[l + n(l - a)]T[l + n(l - 0)]T[1 + n(l - 7)] ' 
that also yields a solution for the electric field in the form of a vector combination 
of F n (x,y, z) functions. 
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There is an essential simplification of the solutions in the planes x = 0, and 
y = 0, and z — 0. One has for these planes 

V- E(a;, i/,*) = -£5—. (42) 
ei - £2 

6. Conclusion 

We demonstrated an application of fractional calculus for electrostatics of 
composite materials in fractal geometry. The method is based on fractional fil- 
tering of Maxwell's equation in the framework of the Gauss theorem by means of 
a convolution integral in Eq. ^ . One should recognize that the coarse graining 
procedure due to the averaging over the period In b is important for the appli- 
cation of fractional calculus. We explore the result of Ref. [20| with subsequent 
development by splitting the composite volume into the fractal volume Vd ~ r d 
with fractal dimension d and the complementary host volume Vh- Integrations 
over these fractal volumes correspond to the convolution integrals and, even- 
tually, to fractional integro-differentiation. Note that this filtering procedure, 
in our case, is the averaging of the Maxwell equation, and not of the electric 
field. The latter is the solution of the obtained averaged equation. It should be 
noted also that the fractal dimension alone is not enough to describe a fractal 
completely. Therefore, although in the analysis above we discuss the Cantor 
sets, the results are valid for any random fractal distribution of charges, since 
we use the fractal dimension only. 

In this connection, we should also admit that embedding the electric charge 
inside any random fractal, as in Eqs. ([5]) and ([7)1. is based on employing the 
characteristic function \ inside integrands. This integration corresponds to the 
convolution integral, as the result of coarse graining. The situation is the same 
for embedding the polarization charge inside a fractal composite. Although, 
instead of the characteristic function, we have its gradient Vx, as shown in 
Appendix B, this leads to the convolution integral (|B. 6|) . as well. 

One should recognize that it is also possible that the permittivity of the 
charged fractal dielectric is less than the permittivity of the host dielectric, 
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e < 0, for example, for biological macromolecules and polymers [27[. In this 
case the argument of the Mittag-Leffler function is positive, and the asymptotic 
behavior of the electric field is absolutely different. One finds easily that the 
limit e — > — docs not correspond to the homogeneous case with e\ = £2- 
Obviously e = is a singular point of Eq. (|23[) . since e is the parameter at the 
higher derivative, and solution ([55)1 is singular at e -> —0 as well. This case 
needs separate consideration. 

We thank Liz Yodim for language editing. This research was supported in 
part by the Israel Science Foundation (ISF) and by the US-Israel Binational 
Science Foundation (BSF). 

Appendix A: Fractional Integration 



Extended reviews of fractional calculus can be found e.g., in [15|, [16|, |2jJ 
Fractional integration of the order of a is defined by the operator 



where a > 0, x > a and the Gamma function T(z) is defined above. Fractional 
derivation was developed as a generalization of integer order derivatives and 
is defined as the inverse operation to the fractional integral. Therefore, the 
fractional derivative is defined as the inverse operator to Q ^"j namely a D™ f(x) = 
a I~ a f(x) and a I" = a D~ a . Its explicit form is 



For arbitrary a > this integral diverges, and as a result of this a rcgularization 
procedure is introduced with two alternative definitions of a D^. For an integer 
n defined as n — 1 < a < n, one obtains the Riemann-Liouville fractional 
derivative of the form 





(A. 1) 




(A. 2) 



jn 

* L D«f(x) = a D a x f{x) = I--f( x ) , 



(A. 3) 



and fractional derivative in the Caputo form 




/(*)• 



(A. 4) 
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There is no constraint on the lower limit a. For example, when a — 0, one has 
^D^x 13 — J r(^+i-a) ^ • This fractional derivation with the fixed low limit is 
also called the left fractional derivative. However, one can introduce the right 
fractional derivative, where the upper limit a is fixed and a > x. For example, 
the right fractional integral is 

= [ a (y ~ ^ 1 fiv)dy ■ (A. 5) 

"(a) J x 

Another important property is 

D a D 13 = D a+I3 and D a I^ = , (A. 6) 

where other indexes are omitted for brevity. We also use here a convolution rule 
for the Laplace transform for < a < 1 

C[IZf{x)] = s- a f(s). (A. 7) 

Note that in physical applications a treatment of the Caputo fractional deriva- 
tive by the Laplace transform is more convenient than the Riemann-Liouville 
one. In the first case, one has to define < n < a integer derivative at the 
initial point, while for the second case one defines n Riemann-Liouville frac- 
tional derivative as initial conditions. Nevertheless, the use of either Caputo or 
Riemann-Liouville fractional derivatives is equivalent for initial value problems 
and leads to the same results [itI, Q] ■ 



Appendix B 

The last term in the l.h.s. of Eq. (1221) corresponds to the induced charge, 
or polarization charge on the fractal, since the electric field is a discontinuous 
function [3J. We also note that, while the notation V r is the divergence in the 
Gauss theorem V r = ^■^r 2 , it is the gradient in the last term V r x( r ) = dX g^ ■ 
Let us consider the TVth step of the fractal construction. It is a union of disjoint 
intervals A^r. In the limiting case one obtains F a — linijv-^oo (J Ajv- Therefore, 
Wj G F a the characteristic function on every closed interval [rj ,rj + Ajv] , is 
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x(Ajv) = 0(r — rj) — Q(r — rj — Ajv)- Differentiation of the characteristic 
function on the intervals yields 

^x(^n) = S(r - r 3 ) ~ 6(r - r 3 - A N ) . (B. 1) 

Therefore, we have for any interval Ajv and at r = rj 

rIr+A N E(r)r 2 V rX (r) = E(r)r 2 - E(r + A N )(r + A N ) 2 . (B. 2) 

This expression is not zero in the limit A at — > 0, since e-zE(r — 0) = E\E(r + 0) 
and £\E{r + A N - 0) = £2E(r + A N + 0) [3]. For brevity's sake the index r 
for the electric field is omitted E r (r) = E(r). We make a shift between points 
rj — and rj + A n — 0. We present the result of integration in Eq. (|B. 2[) in 
the form 

r 2 [E{r - 0) - E(r + A N - 0)] + O(Aiv) = 
r 2 [E(r - 0) - E(r + 0) + E(r + 0) - E(r + A N - 0)] + 0(A N ) . 

Since inside the interval A^r the electric field is continuous, we have E(r + A]y — 
0) « E(r + 0) +i?'(r + 0)Ajv- Now the boundary conditions for the electric field 
at rj can be taken into account: eiEijj — 0) = E\E{r d - + 0). This, finally, yields 

lim r ll E{r)r 2 V rX {r) = £ -±^r 2 E{r) Vr e F a . (B. 3) 

Ajv— s-0 £2 

We have from Eq. (|2U)) that Ql}E r (r)r 2 \J r x{r) is a superposition of random 
values, namely 

llE r {r)r 2 V rX {r) = £ r^fo) . (B. 4) 

This can be approximately considered as an "average" value oIrX( r )Er(r)r 2 — > 
oI"E r (r)r 2 . Indeed, let us present the sum in the form 

£ r 2 E r (r 3 ) = £ f r' 2 E r {r')5{r' ~ Tj ) . (B. 5) 

Hence, one obtains the integration of the electric field with the fractal density 
eF &( r ' ~~ r j)- This corresponds to the fractal volume r a . The next step 

15 



is to consider the i nteg ration in Eq. (|B. 5[) as the convolution integral with the 
kernel (r — r') Q_1 [29[. Eventually, we obtain for the polarization charge term 
in Eq. (|22J 

llE r {r)r 2 VrX{r) -> ^^V r a £r(r)r 2 . (B. 6) 

£2 

Finally, we note that the shift in the intervals A^v is taken between point rj — 
and r j + A jv) ■ One can also make the shift between points rj + and r j + A jy + . 
The difference between electric fields at the points yields [1 — ^]E(rj). Taking 
into account that this difference is opposite to the external normal, we change 
the sign and obtain the same expression as in Eq. (|B. 3ft . 
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